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ABSTRACT 

We use the largest homogeneous sample of globular clusters (GCs), drawn from the ACS Virgo Clus- 
ter Survey (ACSVCS) and ACS Fornax Cluster Survey (ACSE'CS), to investigate the color gradients 
of GC systems in 76 early-type galaxies. We find that most GC systems possess an obvious negative 
gradient in (g-z) color with radius (bluer outwards), which is consistent with previous work. For 
GC systems displaying color bimodality, both metal-rich and metal-poor GC subpopulations present 
shallower but significant color gradients on average, and the mean color gradients of these two sub- 
populations are of roughly equal strength. The ficld-of-vicw of ACS mainly restricts us to measuring 
the inner gradients of the studied GC systems. These gradients, however, can introduce an aperture 
bias when measuring the mean colors of GC subpopulations from relatively narrow central pointings. 
Inferred corrections to previous work imply a reduced significance for the relation between the mean 
color of metal-poor GCs and their host galaxy luminosity. The GC color gradients also show a depen- 
dence with host galaxy mass where the gradients are weakest at the ends of the mass spectrum — in 
massive galaxies and dwarf galaxies — and strongest in galaxies of intermediate mass, around a stellar 
mass of Mi « 10^° Mq. We also measure color gradients for field stars in the host galaxies. We find 
that GC color gradients are systematically steeper than field star color gradients, but the shape of the 
gradient-mass relation is the same for both. If gradients are caused by rapid dissipational collapse 
and weakened by merging, these color gradients support a picture where the inner GC systems of 
most intermediate-mass and massive galaxies formed early and rapidly with the most massive galax- 
ies having experienced greater merging. The lack of strong gradients in the GC systems of dwarfs, 
which probably have not experienced many recent major mergers, suggests that low mass halos were 
inefficient at retaining and mixing metals during the epoch of GC formation. 

Subject headings: galaxies: stellar content - clusters; individual (Virgo, Fornax) - galaxies: elliptical 
and lenticular, cD - galaxies: star clusters - globular clusters: general 



1. INTRODUCTION 

Galactic radial gradients in stellar populations are a 
result of a galaxy's star formation, chemical enrichment, 
and merging histories, and thus can be an important 
discriminant of galaxy formation scenarios. Galaxies 
that form in a strong dissipative collapse are expected 
to have steep gradients in metallicity, as the central re- 
gions retain gas more effectively and form stars more ef- 
ficiently. Thus, in isolation, higher mass galaxies formed 
in this way are expected to have steeper negative metal- 
licity gradients due to their deeper potential wells (e.g., 
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IChiosi fc Carrarol [200l iKawata fc GibsoHl [200l . By 
contrast, in galaxies where merging is a dominant pro- 
cess, radial gradients are expected to weak en due to ra- 
dial mixing that occ i irs during merge rs (iWhitd [19801: 
IBekki fc Shiorol [19991 : lKobavashi|[200l . So if the most 
massive, quiescent g alaxies are the ones mos t shaped by 
major merging (e.g. . [van" der Wei et aD[2009t) . one would 
expect their metal licity gradien t s to be relatively fiat. 
Recently, however, iPipino et al.l ()2010[ ) argue that shal- 
low gradients in massive galaxies can also result from 
lower star formation efficiency and do not necessarily re- 
quire extensive merging. 

The existence of negative optical and near-infrared 
color gradients, where the outer regions are bluer, have 
been well-es t ablished in elliptical and di s k galaxie s (e.g. , 
Franx et all [l989l: iPeletier et all [l990l: iMichardI [20051: 
Wu et al.ll2005tlLiu et al.ll2009[) . and have generally been 
interpret ed as gradients in metallic i ty, or sometimes 
age (e.g.. lKobava"s li fc Arimotol 119991 : iKuntschner et al.l 
[200a [Ra^e et al.ll2008[ ). In lower mass galaxies, how- 
ever, gradients appear to be shallower, nonexistent, or 
even positive. This shows that gradient properties can 
be a function of galaxy mass and perhaps refiects the 
greater diversity in the star formation and evolutionary 
histories of low-mass galaxies. Recent results with large 
samples of galaxies show that while the most massive 
galaxies have shallow or flat color gradients, gradients 
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get increasingly negative toward lower stellar mass until 
M* 3 X IO^^Mq, at which point gradients again be- 
come shallower and even positive (jSpolaor et alj [20091 : 
iTortora et al.l l2010a| ) . For early- type galaxies in par- 
ticular, this has been interpreted as an intrinsic cor- 
relation between gradient and galaxy mass — more neg- 
ative at higher mass — modulated by dry merging at 
higher masses, esp ecially for brightest cluster galaxies 
(|Roche et al.l [2010 ) . 

Nearly all previous studies of stellar population gra- 
dients are of the main stellar body (bulge or disk) of a 
galaxy. Given the complex star formation histories of 
galaxies, the effects of age and metallicity are often dif- 
ficult to disentangle, an d require multiband pho tometry 
and spectroscopy (e.g., iMacArthur et al1l2004[ ). More- 
over, these studies say little about the stellar halo, per- 
haps the oldest galactic component. We thus approach 
the issue of population gradients using a unique tool: 
globular clusters. 

Globular clusters (GCs) are among the oldest stellar 
populations in galaxies, and preserve information from 
the earliest epochs of star formation. Population gradi- 
ents in GC systems have not received much attention, 
but one notable exception was the study of metallicity 
gradie nts in the Milky Way GC system bv lSearle fc ZinnI 
(|1978f) . They showed that although the inner halo GCs 
had a negative gradient, the outer halo GCs had no gra- 
dient, leading them to suggest that the outer halo was 
accreted from dwarf-like fragments. 

In both the Milky Way and nearby galaxies, GCs are 
found to be nearly u niversally old, wit h ages greater 
than 8 Gyr (e.g., [P uzia "ct all I200I Ih ansen et al.l 
[2007t iMarm-Franchet a l. 2009r iWoodlev et al.l 120101 1 
Although in extragalactic systems we are mostly lim- 
ited to broadband colors, the lack of any significant age 
spread in GCs, and the fact that they are generally sim- 
ple stellar populations, allows us to interpret GC colors 
as largely representative of metallicity. 

The color distributions of GCs in massive galaxies are 
often bimodal, and usu ally interpreted as two metallicity 
subpopulatiqns (e.g.. iGebhardt fc Kisslcr-Patig' Il99£ ' 
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a et al.l[200ll : iKundu k WliitmorQ.2001: PcngetaLl 



2006al) (although there is still uncertainty in the trans- 



formation from color to metallicity, see Yoon, Yi & Lee 
2006). Metal-rich (red) GCs are found to have a more 
concentrated spatial distribution than the metal-poor 
(blue) GCs, which results in the total mean color of GCs 
becoming gradually bluer with projected radius (e.g., 
Rhode fc Zepll200lHjOTdan et al.ll2004at iTamura et alj 
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Many studies of massive galaxies have confirmed that 
GC systems taken as a whole have negative color and 
meta l licity gradients (iGeisler et al.lll996HRhode fc Zeofl 
I2OOI [Jordan et all I2004al : ICantiello et al.l 120071) . The 
conventional wisdom, however, has been that indi- 
vidual metal-rich or metal-poor GC subpopulations 
have no color or me tallicity gradients (|Lee et al.l 119981 : 
iRhode fc Zepi [2OOII ). Additional studies of individual 
galaxies, however, have shown that GC subpopulations 
in M49, M87, NGC 1427, and NGC 1399, and nearby 
brightest cluster galaxies do have a slightly negative 
color gradients (IGeisler et al.l 119961 : iForte et al.l 120011: 
iBassino et all 120061: lHarrisll2009allbl V Furthermore, very 
little is known about color gradients in the GC systems of 



dwarf galaxies, whose systems are dominated by metal- 
poor GCs. Similar to population gradient studies of the 
main bodies of galaxies, investigating the color or metal- 
licity gradients of GC systems across a range of galaxy 
mass can provide direct constraints on the formation of 
GC systems and the merging history of their host galax- 
ies. 

In this paper, we present the results from the first 
homogeneous study of color gradients in the GC sys- 
tems of early-type g alaxies. The AC S Virgo Clus- 
ter Survey (ACSVCS. ICote et al.ll200l and ACS For- 
nax Cluster Survey fACSFCS. I Jordan et al.ll2007ah ob- 
served 100 galaxies in the Virgo Cluster and 43 galax- 
ies in the Fornax Cluster using the Hubble Space 
Telescope Advanced Camera for Surveys (HST/ACS). 
All 143 objects are early-type galaxies and range in 
mass from dwarf to giant galaxies. One of the main 
goals of the surveys is the investigation of extragalactic 
GC systems, and pr evious studies have examined their 
color distributions (iPeiig et al.l l2006 al). size distribu- 
tions (iJordan et al.ll2005l: iMasters et a l. 2010 ), luminos- 
ity fu nctions (| Jordan et al.l 120061 I2007bl: IVillegas et al.l 



I2010f ) , formation efficienci es (IPeng et al 
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color- magnitude relations (iMieske et al.l 20061 I2010D . 
Likewise, the surface photometry of the galaxies them- 
selves have also been studied in detail (jFerrarese et al.l 
l2006t iCote et aI]|2007D , allowing us to perform a homoge- 
neous comparison of the color gradients in the field stars 
with those in the GC systems. Another advantage of 
this sample is that distances to most galaxies have been 
determin ed using the method of surface brightn ess fluc- 
tuations ()Mei et al.ll2007t iBlakeslee et al.ll2009l) . Using 
this large and homo genous sample of extragalactic GCs 
(| Jordan et al.l [20091 ) . we measure the color gradients of 
GC systems in the targeted galaxies within the field of 
view (FOV) of the ACS camera, except for four galaxies 
where we use multiple ACS fields. The high resolution 
and quality of the HST images allow us to measure the 
gradients of GC systems in dwarf galaxies as well as in 
individual GC subpopulations for systems showing color 
bimodality. 

This paper is organized as follows: In Section 2, we 
give a description of the GC selection and data analysis. 
The results and discussion are presented in Sections 3 
and 4, respectively. Finally, we conclude in Section 5. 

2. SAMPLE AND DATA ANALYSIS 
2.1. Galaxy Sample and GC Selection 

The data used in this work are drawn from the 
ACSVCS and ACSFCS, which obtained deep, high- 
resolution images for 143 early-type galaxies in the 
F475W (« SDSS g) and F850LP (« SDSS z) filters using 
HST/ACS. These galaxies were selected by morphology 
(E, SO, dE and dSO) and cover a raii ge i n luminosity, 
-22 < Mb < -15 fsee lCote et al.ll2004l and l Jordan et all 
I2007al for details). 

Details about the selection of over 12,000 GC candi- 
dates in the 100 early-type Vi rgo galaxies are described 
in IJordan et all (j2004bl . ]2009l) . Briefly, after selecting 
preliminary GC candidates using magnitude, ellipticity, 
and a broad color cut, all the candidates arc fit with a 
PSF-convolved King model using the KINGPHOT code. 
The probability of an object being a GC (the pcc pa- 
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rametcr) is determined in the plane of magnitude and 
half-light radius with comparison to a number of control 
fields. GCs in the 43 Fornax galaxies were selected us- 
ing the same method. Although previous studies have 
used a criterion of pcc > 0.5 for GCs, in this work, we 
select GCs with pcc > 0.95. The reason we choose this 
stringent criterion is that for the outer regions of dwarf 
galaxies, the contamination from background galaxies is 
the limiting factor. Such a strict selection causes us to 
lose fainter GCs (affecting our completeness), but in- 
creases our efficiency. This stricter cut in the pcc pa- 
rameter actually introduces a varying completeness with 
galaxy mass — essentially, galaxies with more GCs have 
fainter limits — but we have run Monte Carlo simulations 
to show that our results do not change if we choose a 
simple magnitude limit for all galaxies. Our more de- 
tailed but still rigorous approach to selection allows us 
to optimize signal-to-noisc, especially for bimodal color 
distributions where we are splitting the sample in two. 

Contamination by compact background galaxies is 
one of our main problems. To estimate the contam- 
ination from foreground and background, we used 16 
control fields at high latitude fTable 1 of Peng et aD 
2006al). As describe d in detail bv lPeng et all (|2006al ) and 



Jordan et al.l ()2OO90 , the expected contamination was es 



timated for each target galaxy. We have checked that 
the contamination is negligible if we select GCs with GC 
probability pcc > 0.95, averaging ~ 1 object per ACS 
field. 

2.2. GC Subpopulations 

Previous studies show that most massive galax- 
ies have bimodal GC colo r distributions (e.g . 
Gcbhardt & Kisslcr-Patig? 119991: iLarsen et all 120011 : 
Kundu &: Whitmorc. 200l[ ISpitler et all 12005 " 



Peng et all ()2006al ) presented the color distributions 
for GC systems in 100 ACSVCS galaxies. Following 
their work, we use K aye' s Mixture Mod e l (KM M; 
iMcLachlan fc BasfordI 119881: I Ashman et"all Il994f ) to 
decompose the data into two Gaussian distributions 
with the same a (homoscedastic) . We choose the 
homoscedastic case because it is more stable for small 
samples. In practice, for galaxies with large numbers of 
GCs, allowing a to vary has no effect on these results. 
For each GC system, we determine which GCs are 
members of the blue and red GC subpopulations and 
the 'p- value' (not to be confused with pcc) for the 
bimodal model. We consider the GC color distribution 
to be bimodal if the 'p-value' is less than 0.05. A total of 
40 galaxies meet this criterion. Membership in the red 
or blue subpopulation is determined by the membership 
probabilities output by KMM, and corresponds to the 
"dip" in the color distribution. If the p-value is not less 
than 0.05, the galaxy is deemed to have only one GC 
population. 

We show two galaxies as examples in Figure [TJ The 
right panels of Figure [1] show the color distributions for 
GCs in FCC 47 (NGC 1336, panel c) and FCC 153 
(IC 335) panel d). The GC color distribution of FCC 47 
displays two peaks, while the GCs in FCC 153 has just 
one peak in color. The blue and red curves in panel c 
are Gaussian fits to the blue and red GCs determined by 
KMM. In panel d, the black curve shows the best fitting 
of color distribution of whole GC systems using a single 



2.0 



1? 1.5 
I 



1.0 



■FCC0047 Ng,= 208 

:G„|=-0.213±0.033 

■ ^led — U.U / J+U.U3 1 Ot,„,- 


(a) 

— U.U4zj:U.U2j 


■ ■ ■ 1 ■ ■ 


(c) 




^ ^ ' 

v '- '- " ■ 


' ' ■ "Wl 

r 1' 












. ■ ■ ■ ■ • 












r : 





0.1 



1.0 



2.0 



1.5 



I 



1.0 



R (arcmin) 



20 40 60 
N 



■ FCC0153 N^= 34 (b) 

:g,„=-o.ii2±o.o59 ; 


(d) 

> 


i • * '. ' .* ■ ' 





R (arcmin) 



5 10 

N 



Fig. 1. — Color profiles (left panels) and color distributions (right 
panels) of GC systems in two sample galaxies, a) GC system color 
profile of FCC 47 (NGC 1336), a GC system with a bimodal color 
distribution. Each small dot denotes a GC color coded for the blue 
and red subpopulation. The blue, red and black lines arc the best 
linear fit of metal-poor, metal-rich and total GC populations, b) 
GC system color profile of unimodal galaxy FCC 153 (IC 335). 
c) GC color distribution of GCs in FCC 47. Blue and red curves 
represent the Gaussian fitting of blue and red GCs determined by 
KMM. d) GC color distribution for FCC 153, black curve represents 
the Gaussian fitting of all GCs. 

Gaussian function. 

2.3. Calculating Radial Gradients 

The radial gradients are calculated by a linear least- 
squares fit between the GC color or metallicity and the 
logarithm of the radius, defined as: 
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A(g-z) 
Alogi? 
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Alogi? 



(1) 
(2) 



In other words, we measure the change in color or metal- 
licity per dcx in radius. The color gradient errors are 
one-sigma errors and come from linear regression. To 
ensure adequate signal-to-noise, we restrict ourselves to 
the 78 galaxies with at least 20 GCs that meet our se- 
lection criteria (pcc > 0.95). We subsequently eliminate 
VCC 1938 from our sample because of its close projected 
separation from the dwarf elliptical, VCC 1941. We also 
remove the SO/a transition galaxy, FCC167, due to un- 
certainties in measuring its stellar mass. This leaves us 
with a sample of 76 galaxies. If the galaxy has a bimodal 
GC color distribution, we calculate both the color gradi- 
ent of whole GC system and the color gradient of each 
GC subpopulation. We divide the GCs in to red and blue 
using a simple color cut determined by the KMM proba- 
bilities. We tried various approaches, including running 
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KMM as a function of radius, and allowing the dividing 
color to vary with radius using an iterative fitting pro- 
cess. In the former case, the number of GCs limited the 
effectiveness to only a handful of galaxies, and in the lat- 
ter case our results did not change in a significant way so 
we ultimately chose to use the simplest method. Given 
the FOV of ACS (3.'4 x 3.'4), the maximum possible outer 
radius of our gradient measurements is ~ 2.'4, which cor- 
responds to 11.5 and 14.0 kpc at the distances of the 
Virgo and Fornax clusters, respectively. 

The left panels of Figure [T] show color profiles of GC 
systems in FCC 47 (panel a) and FCC 153 (panel b). 
We can see from the figure that both red and blue GC 
systems in FCC 47 show negative color gradients with 
color becoming gradually bluer from the center to the 
outskirts of the galaxy. The gradients of the whole GC 
system is steeper than that of each subpopulation due 
to the dominance of blue GCs at large radii. The GC 
system in the unimodal galaxy FCC 153 also shows a 
shallower but significant negative color gradient. 

For the two most luminous giant galaxies in the Virgo 
cluster: M49 (VCC 1226) and M87 (VCC 1316), the 
gradients are extended by including the GCs in nearby 
ACS fields. Because some targeted galaxies were located 
in the halos of the giants, and th eir own GC syste ms 
appear to be entirely stripped (see iPeng et all 120081 for 
details), we consider the GCs in these fields as part of 
the GC systems of the giant galaxies. For M49, we use 
VCC 1199 and 1192, extending our study to a radius of 
4f5 (22 kpc). For M87, we use VCC 1327, 1297, 1279, 
1185 and 1250, which extends our study to a radius of 
21f3 (102 kpc). There are two similar cases in the For- 
nax cluster. FCC 202 is near FCC 213 (4.'6, 27 kpc) and 
FCC 143 is near FCC 147 (4.'8, 28 kpc). 

3. RESULTS 

3.1. Color Gradients 

Figure [2] presents results for all 76 galaxies, showing 
the color gradient distribution (right panels), and the 
strength of the color gradients as a function of galaxy 
stellar mass (left panels). The ste llar masses for th e 
ACSVCS galaxies were taken from iPeng et al.l (|2008[ ). 
and the masses for the ACSFCS galaxies were calcu- 
lated in the same way as described in that paper us- 
ing g-z photometry from the ACS images (Ferrarese et 
al., in prep) and J- K colors from the Two Micron All 
Sky Survey (;2MASS. ISkrutskie et al.ll2006h . From top 
to bottom, this figure shows the color gradients of red 
GC populations, blue GC and unimodal populations, 
and whole GC systems, respectively. We combine the 
blue GCs and unimodal populations on the same plot 
because unimodal populations for low mass galaxies con- 
sist nearly entirely of blue GCs, and are likely the low 
mass ex tension of the blue GCs in more massive galax- 



studies (e.g. Geisler et al.l 119961: 


Rhode & Zepfl 120011: 


iJordan et al.ll2004al iTamura et al.l 


20061: lHarrifJl2009bn. 



we find that the whole GC systems of most giant early- 
type galaxies have negative color gradients. We also find 
that not only giant galaxies, but also most intermediate- 
and low-mass galaxies show shallow but significant color 
gradients in their GC systems. 
As described in Section 12.21 and 12. 3[ we calculate the 
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Fig. 2. — Color gradients as a function of galaxy stellar mass 
(left panels) and gradient distributions (right panels). From top 
to bottom: red GC systems, blue and unimodal GC systems, and 
all GCs. The dotted lines in panel a, b, and c denote the zero 
gradients. In panel d, e and /, the dot-dashed lines describe the 
mean Gg—z of red, blue and whole GC systems respectively. The 
dashed histogram in panel e is the distribution of color gradients 
in unimodal galaxies. Big filled circles in panel b and c denote the 
mean gradients in given mass bins. The larger and thicker open 
symbols at the high mass ends of panels a b and c, denote the four 
galaxies whose profiles were extended by the use of neighboring 
ACS fields: M49, M87, FCC 213, and FCC 147 (see Section [Z3] l . 



color gradients of individual red and blue GC systems 
respectively if GCs display bimodal color distribution. 
Figure [2] shows that although the subpopulation gradi- 
ents for individual galaxies are often not by themselves 
very significant, both red and blue GC systems have sig- 
nificant shallow negative color gradients when we com- 
bine data from many galaxies. The red and blue GC 
gradients have mean values equaling —0.048 ±0.010 and 
—0.041 ± 0.006 respectively. The errors in the color 
gradients of individual GC systems are taken into ac- 
count when calculating the mean color gradients and 
their errors. Color gradients of red GC populations are 
slightly steeper than that of blue GC systems and seem 
to show more scatter with dispersions a^^d = 0.062 and 
fbiuc = 0.038. Furthermore, both red and blue GC sys- 
tems individually have much shallower color gradients 
than that of whole GC systems (-0.112 ± 0.009) with 
dispersion (Tail = 0.077. 

Table [T] and [5] lists the color gradients of blue, red and 
whole GC systems for our sample galaxies. We only show 
the results for the 76 galaxies with more than 20 GCs 
that meet our selection. 

In panels b and c of Figure [21 big filled circles dis- 
play the mean color gradients in given mass bins with 
bin widths of 0.6 dex. For galaxies with Mi, < 10^°Mq, 
there appears to be a weak correlation between color gra- 
dients and galaxy mass, with color gradients tending to 
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Fig. 3. — The same as panel c and / of Figure [21 but with color 
gradients converted to metallicity for whole GC systems. 

be shallower for dwarf galaxies. But for the higher mass 
galaxies, the trend is flattened, even reversed. In this fig- 
ure, Virgo and Fornax galaxies are plotted together. We 
do not find a significant difference in behavior between 
galaxies in the different clusters. 

3.2. Metallicity Gradients 

Since most GCs are old, single stellar populations, 
trends in their integrated color are generally equated 
with trends in metallicity. Recent studies have found 
a non-linear but mono tonic relation between metallicity 
and color of GCs (e.g., 
I2006at iBlakeslee et al.l 



Harris fc Harri!jl2002l: fPeng et al.l 
20101) . with broadband color less 



sensitive at lower metallicity. IBlakeslee et al.l (j2010D fit 
the color-metallicit y relation from the data shown in 
(jPeng et al.l l2006al) using a quartic polynomial (their 
Equation 1). Although the conversion from color to 
metallicity is still uncer tain and contains con siderable 
scatter, we can use the IBlakeslee et al.l (I2010D relation 
to derive a radial metallicity profile for each GC system. 
After conversion, the metallicity distribution of GC sys- 
tems in many galaxies are not bimodal (see | Yoon et al.l 
[2006t ICantiello fc Blakesleel[2007t IBlakeslee et al.ll2oToir 
but inter preting this is bey ond the scope of this paper 
(see also lSpitler et al.ll2008D . In this work, we only use 
this relation to calculate the mean metallicity gradients 
of the entire GC system of each galaxy. 

Figure [3] presents the distribution of metallicity gradi- 
ents of all GCs and the gradient-mass relation. Sim- 
ilar to the color gradients, the GC systems of most 
galaxies have shallow but significant metallicity gradi- 
ents with a mean value of —0.387 ±0.034 with dispersion 
o'[Fe/H] = 0.284. The metallicity gradient-mass relation 
is also similar to the color gradient-mass relation shown 
in Figure [2] 

3.3. The Gradient-Mass Relation and Gomparison to 
Field Stars 

We have seen that the strength of the GC systems color 
gradient varies as a function of galaxy mass. Similar 
behavior has been seen in the color gradients of the field 
stars for galaxies in other studies. One advantage of the 
ACSVCS and ACSFCS data is that we can also measure 
color gradients for the host galaxies using the exact same 



filters that we use for the GCs, allowing us to make a 
d irect comparison betwe en the two galactic components. 

iFerrarese et al.l ()2006[ ) measured the isophotal light 
profiles of 100 early-type galaxies in the ACSVCS in both 
the g and z band. The light profiles of 43 ACSFCS galax- 
ies are measured by using the same method (Ferrarese 
et al., in prep). The calculations of the color gradients 
of these galaxies are based on their surface photome- 
try. To remove the effect of nucl ei of early-type g alaxies 
(about 2% of effective radius, see lCote et al.li200a [20071) 
and eliminate the significant contamination of sky back- 
ground in outskirts, we measure the color gradients in the 
range from 0.02 Re to Re- The definition of color gra- 
dient is the same as that for GC systems (Equation [Ij . 
We list the steUar color gradients for the ACSVCS/FCS 
galaxies in Table |T] 

Figure m shows the relationship between the field star 
color gradients and galaxy stellar mass, Al*. Larger stars 
denote the mean gradient in a given mass bin. Simi- 
lar to what has been found in other studies, the color 
gradients of the ACSVCS/FCS early-type galaxies are 
mostly negative, with low mass galaxies having flat or 
positive gradients. We find no significant difference be- 
tween color gradients of galaxies in the Virgo and For- 
nax. In the mean, the steepest gradients are found in 
intermedi ate-mass galax i es, wh ich is consistent with the 
finding of iTortora et al.l (l2010al) from SDSS surface pho- 
tometry of galaxies. The circles in Figure U) show the 
mean values for color gradients of GC systems, shown 
in bottom panel of Figure H] We can see that the color 
gradients of GC systems are systematically steeper than 
those of the field stars, but that the gradient -mass re- 
lation is similar in shape to that of stellar systems of 
galaxies. 

At the bright end, the color gradients of GC systems 
seems to be getting shallower again, but there is an im- 
portant caveat. The ACS observations of our sample 
galaxies have field view of 3.'4 x 3.'4. which at the dis- 
tance of the Virgo Cluster is roughly 16 kpc on a side. 
For the massive galaxies, we are only probing the in- 
nermost regions of the halo, ev en when using fields t hat 
observed neighboring galaxies. IRhode fc Zepj (j2001|) . in 
a wide-field study of the GC system in M49 (NGC 4472, 
VCC 1226), found color gradients within 8', but also 
found that th e gradient disa ppeared when expanding the 
radius to 22'. iHarrisI (j2009aD measured color gradients for 
the GCs in M87 and found detectable gradients out to 
~ 60 kpc, or 12f5, a roughly similar radial scale. In this 
work, we have mostly only calculated color gradients of 
the central part of the GC systems of massive galaxies 
due to the limited field of view of ACS, so for massive 
galaxies these color gradients are best described as those 
for "inner halo" GCs. A more comprehensive study of 
GC system color gradients will require wide-field imag- 
ing, such as that being taken for the Next Generation 
Virgo Cluster Survey (Ferrarese et al. 2011, in prep.). 

4. DISCUSSION 

4.1. A Note on Projection Effects 

When we calculate the GC system color gradients, we 
use projected galactocentric distances, not the true three- 
dimensional distances to the galaxy centers. Projecting 
the GC system onto the plane of the sky weakens the 
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TABLE 1 

Color and Metallicity gradients of GC systems of ACSVCS galaxies 



Name Nqc log{Mi,/MQ) Gbiuc G^^d Gaii G[fo/h] Ggai 



(1) 


(2) 


(3) 


(4) 




(5) 


(6) 




(7) 




(8) 


VCC 1226 


950 


11.73 


-0.032 ±0 


017 


-0.037 ±0.019 


-0.139 ±0 


024 


-0.493 ±0 


080 


-0.058 ± 0.003 


VCC 1316 


2002 


11.48 


-0.049 ± 


008 


-0.049 ±0.011 


-0.171 ±0 


013 


-0.570 ±0 


044 


-0.015 ±0.002 


VCC 1978 


709 


11.53 


-0.027 ±0 


027 


-0.021 ±0.023 


-0.070 ±0 


036 


-0.214 ±0 


117 


-0.026 ±0.002 


VCC 881 


256 


11.46 


-0.010 ±0 


036 


-0.027 ±0.044 


-0.011 ±0 


049 


-0.065 ± 


180 


9.999 ± 0.000 


VCC 798 


258 


11.27 


-0.075 ±0 


033 


-0.099 ± 0.033 


-0.150 ±0 


044 


-0.543 ± 


145 


0.104 ±0.006 


VCC 763 


408 


11.37 


-0.048 ± 


025 


-0.145 ±0.034 


-0.117 ±0 


040 


-0.404 ± 


137 


-0.042 ± 0.003 


VCC 731 


772 


11.35 


-0.056 ±0 


024 


-0.033 ±0.021 


-0.118 ±0 


030 


-0.408 ± 


091 


-0.104 ±0.002 


VCC 1535 


163 


10.89 


0.079 ± 


040 


-0.013 ±0.045 


-0.134 ±0 


068 


-0.317 ±0 


227 


0.133 ±0.154 


VCC 1903 


244 


10.92 


-0.065 ±0 


032 


-0.035 ±0.030 


-0.082 ±0 


044 


-0.281 ±0 


132 


-0.102 ±0.004 


VCC 1632 


355 


10.98 


-0.045 ±0 


030 


0.012 ±0.031 


-0.064 ±0 


044 


-0.182 ±0 


128 


-0.056 ±0.004 


VCC 1231 


199 


10.73 


0.006 ± 


030 


0.017 ±0.037 


0.019 ±0 


047 


0.035 ± 


149 


-0.056 ±0.003 


VCC 2095 


75 


10.72 








-0.072 ±0 


067 


-0.256 ±0 


272 


-0.116 ±0.002 


VCC 1154 


132 


10.89 


-0.125 ±0 


045 


0.082 ±0.057 


-0.074 ±0 


067 


-0.363 ± 


218 


-0.019 ±0.002 


VCC 1062 


129 


10.72 


0.013 ±0 


045 


0.032 ±0.053 


-0.130 ±0 


067 


-0.473 ± 


224 


-0.086 ±0.002 


VCC 2092 


52 


10.68 


-0.093 ±0 


067 


0.095 ±0.067 


0.097 ±0 


103 


0.163 ±0 


345 


-0.132 ±0.003 


VCC 369 


128 


10.51 


-0.002 ±0 


048 


0.027 ±0.055 


0.073 ±0 


075 


0.170 ±0 


230 


-0.078 ±0.003 


VCC 759 


112 


10.65 


-0.056 ±0 


035 


0.023 ± 0.042 


-0.048 ± 


052 


-0.221 ±0 


187 


-0.060 ±0.006 


VCC 1692 


93 


10.53 


0.011 ±0 


035 


-0.073 ±0.057 


-0.094 ±0 


077 


-0.205 ±0 


267 


-0.077 ±0.006 


VCC 1030 


118 


10.12 


-0.075 ±0 


037 


-0.094 ± 0.048 


-0.209 ±0 


064 


-0.669 ± 


217 


-0.081 ±0.009 


VCC 2000 


148 


10.38 


-0.057 ±0 


030 


-0.072 ±0.074 


-0.188 ±0 


045 


-0.684 ± 


171 


-0.092 ±0.002 


VCC 685 


125 


10.49 


-0.035 ±0 


042 


-0.087 ±0.052 


-0.210 ±0 


058 


-0.734 ± 


208 


-0.026 ±0.008 


VCC 1664 


104 


10.42 








-0.095 ±0 


062 


-0.328 ± 


175 


-0.118 ±0.003 


VCC 654 


23 


10.36 








-0.036 ± 


121 


-0.090 ± 


600 


-0.056 ±0.005 


VCC 944 


62 


10.46 


-0.060 ± 


055 


-0.023 ±0.061 


-0.162 ±0 


088 


-0.560 ± 


319 


-0.075 ±0.002 


VCC 1720 


42 


10.31 


-0.154 ±0 


060 


-0.149 ±0.069 


-0.178 ±0 


105 


-0.733 ± 


388 


-0.118 ±0.002 


VCC 355 


29 


10.20 








-0.345 ± 


131 


-1.455 ±0 


586 


-0.083 ±0.005 


VCC 1619 


44 


10.24 








-0.075 ±0 


091 


-0.501 ±0 


367 


-0.026 ±0.001 


VCC 1883 


43 


10.22 


-0.079 ±0 


046 


-0.021 ±0.080 


-0.144 ±0 


067 


-0.452 ± 


207 


-0.018 ±0.004 


VCC 1242 


78 


10.18 


-0.029 ± 


036 


-0.117 ±0.040 


-0.147 ±0 


057 


-0.490 ± 


177 


-0.062 ±0.002 


VCC 784 


43 


10.23 








-0.001 ±0 


114 


-0.082 ± 


406 


-0.074 ±0.003 


VCC 1537 


25 


10.01 








-0.342 ± 


065 


-1.068 ±0 


294 


-0.058 ±0.003 


VCC 778 


43 


10.26 








-0.137 ±0 


099 


-0.387 ±0 


382 


-0.064 ±0.003 


VCC 1321 


22 


9.84 








-0.085 ± 


078 


-0.630 ± 


426 


-0.084 ± 0.004 


VCC 828 


48 


10.14 


0.032 ± 


048 


0.065 ±0.107 


-0.159 ±0 


077 


-0.553 ±0 


365 


-0.038 ±0.002 


VCC 1630 


29 


10.06 








-0.199 ±0 


118 


-1.005 ±0 


479 


-0.060 ±0.002 


VCC 1146 


53 


9.94 








-0.139 ±0 


062 


-0.105 ±0 


105 


-0.081 ±0.003 


VCC 1025 


58 


10.33 








-0.210 ±0 


055 


-0.801 ±0 


239 


-0.113 ±0.003 


VCC 1303 


37 


10.02 








-0.038 ± 


042 


-0.271 ±0 


262 


-0.113 ±0.003 


VCC 1913 


36 


10.03 








-0.037 ±0 


070 


-0.235 ±0 


275 


-0.079 ± 0.004 


VCC 1125 


39 


9.91 








-0.080 ± 


064 


-0.382 ± 


302 


0.042 ± 0.004 


VCC 1475 


52 


9.89 








-0.013 ±0 


059 


0.067 ± 


266 


-0.044 ± 0.004 


VCC 1178 


58 


9.85 








-0.102 ±0 


063 


-0.450 ± 


228 


-0.005 ± 0.004 


VCC 1283 


36 


9.96 








-0.113 ±0 


099 


-0.533 ± 


324 


-0.047 ±0.003 


VCC 1261 


22 


9.69 








-0.033 ± 


070 


-0.202 ±0 


235 


-0.007 ±0.004 


VCC 698 


83 


9.98 








-0.068 ± 


051 


-0.294 ± 


235 


0.018 ±0.006 


VCC 1910 


34 


9.32 








-0.158 ±0 


074 


-0.703 ± 


315 


-0.017 ±0.005 


VCC 856 


24 


9.35 








-0.038 ± 


065 


-0.173 ±0 


274 


0.001 ±0.006 


VCC 1087 


43 


9.52 








-0.110 ±0 


045 


-0.383 ± 


228 


-0.057 ±0.010 


VCC 1861 


28 


9.46 








-0.253 ±0 


107 


-0.974 ± 


479 


-0.009 ±0.007 


VCC 1431 


40 


9.34 








-0.078 ±0 


054 


-0.316 ±0 


196 


0.065 ±0.012 


VCC 1528 


28 


9.21 








-0.124 ±0 


055 


-0.616 ±0 


333 


-0.046 ± 0.006 


VCC 2019 


20 


9.01 








-0.021 ±0 


062 


-0.364 ± 


380 


-0.084 ±0.008 


VCC 1545 


27 


9.15 








-0.122 ±0 


053 


-0.341 ±0 


266 


-0.143 ±0.005 


VCC 1407 


22 


9.09 








0.040 ± 


050 


0.178 ±0 


207 


0.011 ±0.004 


VCC 1539 


24 


8.72 








-0.079 ± 


073 


-0.363 ± 


359 


0.113 ±0.019 



The name of galaxies 
2 Total number of GCs those p > 0.95 

Logarithm of stellar mass (in unit of Mq) 
^ Color gradient of red GCs with error, if bimodal 
^ Color gradient of blue GCs with error, if bimodal 
^ Color gradient of all GCs with error 
^ Metallicity gradient of all GCs with error 
^ Color gradient of galaxy with error 

Note ^Yie FOV of ACS, the outer boundaries of our measurements of gradients in most galaxies are about 2—3 arcmin. 
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TABLE 2 

Color and Metallicity gradients of GC systems of ACSFCS galaxies 



Name A^gc log(A1*/X0) Gbiuc Crcd Gaii Gifo/h] Ggai 



(1) 


(2) 


(3) 


(4) 




(5) 




(6) 




(7) 




(8) 


FCC 21 




1 ± .oo 






n 994 -1- 


041 


n 941 -1- 


u^y 


n 7Qfi -u n 




074 -1- 01 






11 42 




m Q 


041 -1- 


m Q 


n oRfi -1- 




n '^04 -1- 


070 


n ni -1- 009 


FCC 91 Q 


297 


11.14 


— n 047 -1- 


038 


—0.119 lb 


027 




048 


_n fici _|_ f) 

— U.UO-L ^ U 


154 


_n ni q _|_ n n09 
— u.U-Ly ^ U.UU.Z 


Nnc 1 '?4n 

i> VJv_y ± •J'-tKJ 


151 


10.99 


f) nno _|_ n 
w.vjvjo _i_ yj 


037 


—0.044 lb 


066 


— f) 094 H- 


048 


— D 080 H- 

\J.\JU\J _1_ \J 


219 


—0 07*^ + 00'^ 


VCC 97fi 


280 


10.67 


— n o'^n -1- 

— u.uou ^ u 


027 


—0 09Q -1- 

— U.U.ZC/ ^ U 


032 


—0.176 lb 


037 


—0 fi09 -1- 


118 


—0 0fi7 -1- 00"^ 


FCC 147 


264 


10.69 


— 047 -1- 


022 


n noe; _|_ n 


039 


—0.154 lb 


029 


—0 "^1 -1- 


090 


—0 049 -1- 00"? 




Q7 




n 097 -1- 


OAR 


n 1 7*3 _|_ n 


uuu 


n 1 71 -1-0 


074 


n c;i -1- 


232 


n 1 _|_ n 004 

— U.IZO ^ U.UU'i 


FCC 83 


217 


10.51 


— 009 -1- 


021 


_n 1 no _|_ n 


031 


— ORQ -1- 


037 


—0 94J^ -1- 

— U.Z10 ^ u 


115 


—0 1 O'^ -1- 009 
— u.±uo ^ u.uuz 


FCC 184 


230 


10.82 


-0.037 ±0 


046 


-0.089 ± 


048 


-0.013 ±0 


079 


-0.074 ± 


248 


-0.011 ±0.001 


FCC 63 


163 


10.43 


-0.045 ±0 


035 


-0.082 ± 


040 


-0.171 ±0 


046 


-0.513 ±0 


167 


-0.120 ±0.006 


FCC 193 


25 


10.48 










0.023 ± 


140 


0.056 ± 


475 


-0.107 ±0.003 


FCC 170 


44 


10.29 










-0.137 ±0 


074 


-0.578 ± 


350 


-0.007 ±0.005 


FCC 153 


33 


9.94 










-0.112 ±0 


059 


-0.502 ±0 


340 


0.078 ±0.012 


FCC 177 


45 


9.82 










-0.097 ±0 


050 


-0.401 ± 


291 


0.200 ±0.014 


FCC 47 


206 


9.97 


-0.043 ±0 


021 


-0.074 ± 


031 


-0.208 ±0 


033 


-0.687 ± 


116 


-0.093 ± 0.004 


FCC 190 


106 


9.87 


-0.087 ±0 


022 


0.023 ± 


064 


-0.116 ±0 


030 


-0.619 ±0 


156 


-0.005 ± 0.004 


FCC 249 


115 


9.99 










-0.115 ±0 


041 


-0.446 ± 


165 


-0.083 ±0.007 


FCC 148 


58 


10.03 


0.006 ± 


043 


0.019 ±0 


049 


0.042 ± 


076 


0.161 ±0 


305 


0.189 ±0.010 


FCC 255 


53 


9.56 


-0.035 ±0 


034 


-0.044 ± 


041 


-0.147 ±0 


044 


-0.526 ±0 


179 


0.012 ±0.005 


FCC 277 


22 


9.88 










-0.300 ± 


119 


-0.995 ±0 


501 


-0.055 ±0.003 


FCC 182 


30 


9.40 










0.080 ± 


078 


0.204 ± 


278 


-0.036 ±0.010 



The name of galaxies 
Total number of GCs those p > 0.95 
Logarithm of stellar mass (in unit of A4q) 
Color gradient of red GCs with error, if bimodal 
Color gradient of blue GCs with error, if bimodal 
Color gradient of all GCs with error 
Metallicity gradient of all GCs with error 
^ Color gradient of galaxy with error 

Note .(.Q ^Yie FOV of ACS, the outer boundaries of our measurements of gradients in most galaxies are about 2-3 arcmin. 




0.1 - 



0.0 -■ 



-0.1 - 



<f gradient of individual Virgo galaxy 
.Q2 -'^ gradient of individual Fornax galaxy 
^mean gradient of galaxies 
# mean gradient of GC systems 
I I I ■ ■ 



10 

log (M,.,u^,„ 



11 



12 



Fig. 4. — Color gradients of galaxies and GC systems binned by mass as a function of stellar mass, A^*. Small open stars with error bars 
are the gradients and errors in individual galaxies. Filled stars and circles describe the mean color gradients and errors of galaxies and GC 
systems in given mass bins. 
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Fig. 5. — A simulated GC system color profile. Solid line is the 
intrinsic profile, dashed and dotted lines denote color profiles of 
red and blue GC populations when projected onto the plane of the 
sky. 

measured gradients, as GCs projected onto the center 
are actually a mix of GCs at all radii. This is less of 
a problem for more centrally concentrated systems (i.e., 
more steeply rising density profiles toward the center). 
Because the distribution of red GCs can be more con- 
centrated than blue GCs, the projection effects for the 
two subpopulations could be different. 

In order to test the affects of projecting real gradi- 
ents into the plaii e of the sky, we provide one test case. 
ICote et all (I2001D deprojected the spatial distribution 
of the red and blue GC subpopulations in the galaxy 
M87. The projections of the model density distribution 
are consistent with the observed surface density of GCs. 
They obtained: 

= ( 2di^)"(^ + 2oii^)"' ^ < ^^^'^ 

We assume that the initial color gradient of GC systems 
are —0.15 and project the model density distribution into 
two dimensions. The resulting projected color profiles 
are shown in Figure [3 The projection effect flattens the 
gradient of both red and blue GC systems to —0.134 and 
—0.132, respectively. Because the blue GCs are more ex- 
tended than the red ones, the flattening is slightly more 
obvious in the gradient of blue GC system. The total ef- 
fect of projection in this case, however, is relatively small, 
roughly 12%, and the relative effect between the red and 
blue GCs is negligible. We do not correct for projection 
effects because we do not know the three-dimensional 
density profiles of the GC systems. We simply note that 
the true radial color gradients for these GC systems are 
slightly steeper than the projected gradients, but that for 
a realistic density profile of GCs, this correction is likely 
to be of the order of ~ 10%. 

4.2. Gradient- Induced Bias in the Colors of GC 
Subpopulations 



Our results show that the GC systems of early- 
type galaxies display significant negative color gra- 
dients, which is c o nsistent with previ o us work (e.g., 
Strom et al.l [Tosl IGeisler et al.l [l99fil: IRhode fc Zepfl 
20041 : [Jordan et al.l [200'4al : iCantiello et al.ll2007f ). How- 
ever, there exists more uncertainty about whether in- 
dividual GC subpopulations display color gradients or 
not. Some previous studies have found shallow gradi- 
ents in individual G C subpopulations (e.g.. lBassino et ahl 
200fiHH arris"2009a','l?) while other studies hav e not (e.g., 
Forbes et al . 2004; Canticllo ct al. 2007; Ku ndu fc Zepj 
2007f ). In this work, we find that only a few of the in- 
dividual GC subpopulations show significant gradients 
(> 3(t), but the overall trend is obvious when measured 
over 39 galaxies, i.e., individual GC subpopulations dis- 
play negative gradients statistically. It is the first ho- 
mogeneous study of color gradients of GC systems in 
early-type galaxies, and emphasizes the power of using a 
large, homogeneous sample of galaxies. 

One consequence of this result could be on the mean 
measured colors of the subpopulations. The studies with 
the highest precision photometry and leas t contamina- 
tion have often used H ST imaging (e.g., iLarsen et al.l 
120011 : IPeng et aLl l2006aD which necessarily has a small 
field of view relative to the largest nearby GC systems. 
Results from these studies have shown that the mean 
color of the blue and red GC subpopulations is a function 
of galaxy mass or luminosity, where more massive galax- 
ies have redder GCs, in the mean. This has generally 
been interpreted as a mass-metallicity relation for GC 
systems (as opposed to for individual GCs). The corre- 
lation for metal-poor GCs, although weaker than that for 
the red GCs, has drawn particular interest because it im- 
plies a connection between the earliest forming GCs and 
the final halos in whic h they reside (iLarsen et al, 2Q0it 
Burgarella et al.l 120011: IStrader et a"Ln2004l: iPeng et al I 
2006al) . 

Given the fixed and relatively small field of view for 
the instruments, there is a significant aperture sam- 
pling e ffect that varies ac ross the studied range of galaxy 
mass. IPeng et aLl (|2008f ) showed that for galaxies with 
Mb > —18, the entire GC system typically fits within the 
ACS field of view. At higher luminosities, the ACS field 
will miss some fraction of the outer regions. This frac- 
tion can be as high as ^ 90%, in the case of M87. This 
bias toward the centers of galaxies would not matter if 
the GC subpopulations did not possess color gradients. 
We have found, however, that they typically have gra- 
dients of 0.04-0.05 mag/dex in {g-z). This results in a 
bias where the most massive galaxies are sampled where 
the GCs are most red. This would particularly affect 
the blue GCs, which can have a more extended spatial 
distribution. 

Although the resolution of this issue will ultimately re- 
quire precision wide-field imaging, we can estimate the 
degree of bias that color gradients may have introduced. 
For the most massive galaxies in Virgo, such as M49, 
M87, and M60, the effective radii (i?e,GCs) of the GC 
populations are 42, 41, and 24 kpc, respectively. This 
was determined using Scrsic fits to the GC spatial dis- 
tribution from ACSVCS data and publi shed photometry 
(IMcLaughlinlll999HRhode fc Zepjl200l]) . The mean pro- 
jected radius for the GCs observed in the ACS field, and 
for which the subpopulations colors were measured in 
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Fig. 6. — Mean colors of red (circles) and blue GCs (diamonds) 
as a function of galaxy luminosity (Mb) with data from Peng et al. 
(2006) (filled points) and inferred corrections to a common radius 
of -Re^GCs (open points). The correction is only important for the 
most luminous galaxies. The dotted lines show linear fits to the 
data with inferred correction. The slope of the blue GCs is -0.007, 
only half the value when fitting the original data. 



iPeng et all (|2006a[) . is roughly 5 kpc. Given the gra- 
dients for the red and blue populations measured in 
this paper (Tabled]), we can infer the expected differ- 
ence between the mean color in the ACS field and the 
mean color at li?e,GCs, which should roughly represent 
the mean color of the entire GC subpopulation if the 
color gradient is constant at all radii. The color differ- 
ence, Ag^,{ACS - Re,GCs), for M49, M87, and M60 is: 
0.031 mag, 0.046 mag, and 0.019 mag, respectively. 

Could such a shift to the blue affect the previously 
published correlations between galaxy mass and GG sub- 
population metallicity? We estimate Ag-z as a function 
of galaxy Mb in the A CSVCS sample for e asy compari- 
son with the analysis in iPeng et al.l ()2006a|) . We use the 
measured mean colors within the HS T/AGS field of view 
from Table 4 of IPeng et all (l2Q06aD . nc,^ for th e GC 
systems in the ACSVCS galaxies (jPeng et al.ll2Q08l: Peng 
et al., in prep), and the mean color gradients for the 
red and blue GC subpopulations (—0.048 and —0.041, 
respectively. Figure [2]) to infer the mean color for the red 
and blue subpopulations at R^^gcs, which we take to be 
representative of the entire population. We weight each 
galaxy's c ontribution to An _z by their total number of 
GCs from lPeng et all (120081) . 

Figure m mirrors Figure 8a of IPeng et all (|2006aD and 
plots the mean colors versus galaxy M^, both measured 
with AGS (solid points) and inferred at -Re.GCs (open 
points). As expected, the correction is only important 
for the two most luminous bins in Mb- Whereas the 
relation for the red GCs is not significantly affected, as 
it was originally quite steep, the slope for the blue GC 
relation is noticeably smaller. The fit to the original 
measurements gave a blue GC slope of —0.0126 ±0.0025, 



a nonzero slope at the 5cr level. The fit to the newly 
inferred mean colors produces a slope of — 0.0069±0.0025 
(systematic errors from the correction process arc not 
included). This shallower slope is now only significant 
at 2.8ct, and could potentially be even less significant if 
more accurate color gradients on the outer regions are 
measured to be steeper than what we have measured. 

We want to emphasize that this exercise is far from 
conclusive, and only serves to warn that color gradients 
in the GC subpopulations will potentially affect conclu- 
sions drawn from imaging the central regions of galaxies. 
Similar biases have been noted for the color-magnitude 
relation of early-ty pe galaxies wher e colors are measured 
in fixed apertures (|Scodeggioll200lD . The correction that 
we infer out to i?e,GCs assumes a constant color gradi- 
ent over the entire GG system. This assumption is un- 
verified, and probably provides an uppe r limit on the 
possib le co rrection given the results of iRhode fc Zepj 
(|200H) andlHariid (|2009bf ) who find a flattening gradient 
at large radius. Nevertheless, a shallower (or potentially 
non-existent) relation between the mean color (metallic- 
ity) of metal-poor GCs and their host galaxy mass will 
have implications for understanding the formation of GC 
systems, and the solution t o this problem await s high- 
quality wide-field data (e.g.. IRhode fc ZeDj[200l) . 

4.3. The Formation of GC systems and Their Hosts 

That GG systems should have negative color (metal- 
licity) gradients is perhaps not surprising given that GC 
formation is by its very nature a product of high star 
formation efficiency. Most models predict that high 
efficiency of star formation plus metal retention leads 
to more enriched populations at the centers of galax- 
ies. Interestingly, even though GCs are among the old- 
est objects in galaxies, and thus have presumably ex- 
perienced the largest amount of merger-induced radial 
mixing of any stellar population, the color gradients in 
most intermediate- and high-mass galaxies arc still sig- 
nificantly negative. 

Di Mattco e t al.l ()2009f ) investigated the survival of 
metallicity gradient after a major dry merger. For ellip- 
ticals with similar initial gradients, they concluded that 
the final gradient is about 0.6 times of the initial after a 
major dry merger. Dissipational mergers, however, can 
either fiatten or enhance gradients, depending on the ini- 
tial gradients and the amount of gas involved. This de- 
pendence on merger history is one of the reasons why 
gradients in massive galaxies are expected to have larger 
dispersion. We have very few galaxies on the high-mass 
end, so it is difficult for us to probe the dispersion in GC 
system color gradients in this mass range. For the high- 
mass galaxies in our sample, we are also only probing 
the very inner halo, so it is possible that the gradients in 
this region arc either more robust to dilution or had the 
strongest initial gradients. The color gradients for blue 
GCs (presumably the oldest GCs) in massive galaxies are 
detectable but shallower than at lower masses, and this 
may be a sign that mergers have played a part in their 
history. It would be interesting to extend this study to 
wider fields of view for the more massive galaxies. 

We find that the mean gradients for the red and 
blue subpopulations arc similar in magnitude (—0.048 
and —0.041, respectively), but its interpretation is 
confounded by the uncertain conversion from (5 — 
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z) to metallici ty. Both the iPeng et ahl (|2006aD and 

iBlakeslee et aLl (|2010[ ) transformations have slopes that 
are roughly 3-4 times steeper at the mean blue GC color 
than at the mean red GC color. By extension, the true 
metallicity gradient for metal-poor GCs should be 3-4 
times steeper than that for metal-rich GCs (given their 
similar gradients in color). This would be a fairly re- 
markable difference between the two populations, but is 
still entirely dependent on the assumed color-metallicity 
relation. We plan to revisit this question when the trans- 
formation from (g-z) is better understood. 

The relationship between the color gradients of GC 
systems and host galaxy mass offers some interesting in- 
sights into the formation and evolution of stellar halos 
in early- type galaxies. Even for measurements of color 
gradients from galaxy surface photometry, it was only 
relatively recently that the data quality and mass range 
prob ed has been su fficient t o study trends in galax y mass 
(e.g. iForbes et al.ll2 005: Sp olaor et aL|[200 9: Rawl c et al.l 
[20ia iTortora et al.li2010ai) . iTortoraetaLi ()2010bD have 
also run simulations to show that environment can also 
play a role in the observed gradients. Our results show 
that the shape of gradient-mass relation for GC sys- 
tems is similar to that for the galaxies themselves, with 
a minimum around ~ IO^^Mq (Figure[5]). That the GC 
color gradients are universally steeper than those for the 
field stars is an interesting result. If the GCs formed 
in higher efficiency s tar forming event s than the bulk of 
the field stars (e.g., iPeng et all 120081 ). then that might 
result in steeper gradients. One caveat for the inter- 
pretation, however, is that the total GC gradients are 
actually a combination of the red and blue GC popula- 
tions, which may not have formed contemporaneously in 
the present-day halo. The steeper gradients are likely a 
combination of the increasing fraction of blue GCs and 
the in creasing specific freque ncy of GCs at lower metal- 
licity (jHarris fc Harris! |2002|) . The color gradients for 
the individual GC populations are similar in magnitude, 
if not slightly shallower than the gradients for the field 
stars. 

The shape of the gradient-mass relation for both GC 
systems and field stars is broadly consistent with a 
model where color (metallicity) gradients are increas- 
ingly steeper in higher mass halos due to metal reten- 
tion, but then are diluted in the highest-mass galax- 
ies {Mi, > 10^°7V(o) due to the increasing importance 
of mergers in their evolution. One difference between 
the GCs and the stars is that the stars in some dwarfs 
exhibit significantly positive color gradients, which are 
often interpreted as d ue to age gradients (age iiicreas- 
ing with ra dius) fc.g.. iLa Barbera fc de Carvalhol [20091 : 
iSpolaor et al. 2010,). This is not difficult to produce if 
there has been recent low level star formation at the 
galaxy center. However, we do not see any case of this 
for the GC systems, nor might we expect to as the star 
formation rate density required to produce young GCs is 
much higher than needed to produce a slight age gradient 
in the field. We notice that there is a prominent outlier 
in Figure H VCC 798 (M85) with mass - IO^ ^^TWq, 
which has a steep positive color gradient. This galaxy 
is known as a young, gas-rich merger rem nant (e.g., 
iSchweizer fc Seitzedll99alPeng et al.li2006bl) and hosts a 
large-scale stellar disk ( Ferrarese et al.ll2006f ) . During the 
gas-rich merger, the central starburst produced a young. 



blue stellar population in the center of galaxy. There- 
fore, the positive color gradients are cominon in gas-rich 
merger remnant (e.g.. lYamauchi fc Gotol[2005D . How- 
ever, the color gradient of the GC system in VCC 798 is 
negative and quite normal. One of the possible reasons 
is that the number of GCs formed during the merger is 
negligible compared to the preexisting old GC popula- 
tion. The fact that the GC systems of dwarf galaxies 
have shallow or flat color gradients suggests that metal 
retention and mixing were not efficient during the epoch 
of GC formation. 

5. CONCLUSION 

We use HST imaging from the ACS Virgo and For- 
nax Cluster Surveys to conduct the first large-scale study 
of globular cluster system radial color gradients. We 
present results for 76 early-type galaxies, measuring [g- 
z) color gradients for GC systems across a range in galaxy 
stellar mass (8.7 < log(X*/X0) < 11.8). For 39 galax- 
ies whose GC systems show significantly bimodal color 
distributions, we also measure the color gradients in the 
GC subpopulations. Using the surfac e phot ometry of 
ACSVCS galaxies from iFerrarese et al.l (|2006D . we mea- 
sure the radial color gradients of the field stars in the 
same galaxies and same filters, allowing a direct compar- 
ison of GC and field star radial gradients. Wc caution 
that the FOV of ACS means we only measure the cen- 
tral part of many large galaxies, which may introduce 
an aperture bias if the color gradient of galaxies are not 
constant with the radius. We find that: 

1. GC systems as a whole have negative color gradi- 
ents, with an average gradient over the entire sam- 
ple of —0.112 ± 0.009 mag in {g — z) per dcx in 
radius. 

2. On average, red and blue GC subpopulations also 
show significantly negative color gradients at the 
mean level of -0.048 ± 0.010 and -0.041 ± 0.006, 
respectively. Although a gradient is sometimes 
difficult to detect for any individual galaxy, the 
combined sample shows this property with higher 
signal-to-noise. 

3. We find a relationship between GC system gradient 
strength and galaxy stellar mass, where the gradi- 
ents are fiat at low mass, increasingly negative with 
mass until M.^, « IQ^^Mq and then staying con- 
stant or less negative at higher mass. This trend 
parallels the gradient-mass relationship we find for 
the field stars in the ACSVCS galaxies. The GC 
system gradients are systematically steeper than 
that for the field stars, which is likely a reflection of 
the dominance of blue GCs at large radius. These 
observed trends, however, are limited by the small 
number of galaxies at high and low mass in our 
sample. 

4. Color gradients in the GC subpopulations can 
cause a bias in the measurement of the mean colors 
of GCs when the data only covers the central region 
of the galaxy. We infer a correction using the mea- 
sured gradients and find that the slope between the 
mean color of metal-poor GCs and the luminosity 
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of their hosts can be reduced by nearly a factor of 
two from previous measurements, raising questions 
about its level of significance. 

5. The shape of the gradient-mass relation for GC 
systems is consistent with picture where the for- 
mation and chemical enrichment of the GC sys- 
tem becomes more efficient as the mass of the host 
galaxy increases, but is further affected by signifi- 
cant merging and radial mixing in the most massive 
galaxies. 

6. In a test case, the intrinsic, three-dimensional color 
gradients are likely to be roughly 10% steeper 
than the projected gradients given a reasonable 
spatial distribution of GCs. 
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